Approximate analytical results on the cavity dynamical Casimir effect in the presence 

of a two-level atom 
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We study analytically the photon generation from vacuum due to the Dynamical Casimir effect 
in a cavity with a two-level atom, prepared initially in an arbitrary pure state. Performing small 
unitary transformations we obtain closed analytical expressions for the probability amplitudes and 
other important quantities in the resonant /dispersive regimes. 
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I. INTRODUCTION 

The fascinating dynamical Casimir effect (DCE), i.e., 
the creation of quanta from vacuum due to the motion of 
macroscopic neutral boundaries (or due to the time vari- 
ations of material properties of these boundaries, such as 
the dielectric constant or conductivity), attracted atten- 
tion of many theoreticians for several decades since the 
first publications [l|, Q (see for the most recent re- 
views) . Quite recently the first experiments on modelling 
this effect in the superconducting stripline waveguide ter- 
minated by a SQUID subjected to rapidly varying mag- 
netic flux (resulting in time-dependent boundary condi- 
tions simulating the motion of some effective boundary) 
were performed This realization can be called the 
"single mirror DCE" 0, 0| ■ One of its specific features 
(which was used as one of decisive proofs of the effect) is 
the creation of correlated photon pairs emitted outside. 

Another possible realization corresponds to the case of 
a closed cavity with moving wall(s). This "cavity DCE" 
was considered for the first time in [l|, and it attracted 
the special attention, because the number of photons ac- 
cumulated inside the cavity can be significantly increased 
in the case of periodical motion of the wall(s) under cer- 
tain resonance conditions fi-[i3]. The simplest Hamiho- 
nian describing this effect in the absense of dissipation 
reads (we assume h = 1) 9] 
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uit)n + ixma^^ ~a^), x{t) ^ 



duj/dt 



(1) 



where a and are the cavity annihilation and creation 



operators, n 



'a is the photon number operator, and 



iij{t) is the cavity "instantaneous" eigenfrequency, which 
depends on time due to the time-dependent geometry 
of the cavity. This Hamiltonian transforms the initial 
vacuum state into the squeezed vacuum state, so that only 
even numbers of quanta can be generated with nonzero 
probabilities. Several possible realizations of the cavity 



DCE Hamiltonian ([T]) were proposed a few years ago 111, 
[l^ . and the experimental progress was reported in 13 1 
(other schemes, based on the fast optical modulation of 
the cavity length, were proposed recently in In 
view of this progress, the problem of detecting the created 
photons becomes a timely one. 



One of the simplest detectors could be a two-level sys- 
tem ("atom") [il, which can serve as an approxi- 
mate model of either real Rydberg atoms (or bunches of 
such atoms) [l^, 113 oi' some kinds of "artificial atoms" 
(made, e.g., from Josephson's contacts used in quantum 
superconducting circuits |18| - in this case one deals 
with the "circuit DCE" [3, [T^). In such a case, one 
should add to Hamiltonian ([1]) the free atom Hamiltonian 
Ha = (fl/2)(Tz and the interaction Rabi Hamiltonian 



Hp 



g {a + a^) (cr+ + a-) 



(2) 



where a± and az are the standard Pauli operators. 



|e)(e|-|g)(5l, 



l5)(e|, a+ = \e){g\. 



il and g are the atomic transition frequency and the 
atom-field coupling constant (assumed real), respectively. 
The kets \g) and |e) can be interpreted as "atomic" 
ground and excited states, respectively. 

There are two possible scenarios. In the first one the 
atoms are injected into the cavity after the walls made 
a sufficient number of oscillations and returned to initial 
positions. Then the solution is splitted in two steps: first 
one calculates how many photons are created in the cav- 
ity, using Hamiltonian ([1]), and after that one turns on 
the interaction between the atoms and the field, using 
Hamiltonian ([2]) together with the free field Hamiltonian 
with UJ — const. This case was analyzed recently in 

Here we study another case: when the detecting atom 
is present in the cavity for all time, infiuencing the pho- 
ton generation process. This situation looks quite realis- 
tic, because the cavities in the experiments proposed in 
[ll|,[l3] have the dimensions of the order of a few centime- 
ters (with the resonance frequencies of the order of a few 
GHz, corresponding to the transitions between nearest 
levels of highly excited Rydberg atoms). The travel time 
of atoms with velocities of the order of 100 m/s through 
such cavities is of the order of 10^"' s (or bigger for slower 
atoms). On the other hand, the time of oscillations of 
the cavity wall (more precisely, some effective "moving 
plasma mirror" created by periodic laser pulses illumi- 
nating a thin semiconductor slab attached to the wall) 
in the experiments discussed is expected to be of the or- 
der of 10^^-10^^ s, so that it is much shorter than the 
atom travel time. Therefore it could be more easy to 



2 



send atoms through the cavity continuously, instead of 
adjusting the exact moment of their entrance into the 
cavity. In such a case, one can assume that the atom is 
permanently inside the cavity, thus interacting with the 
field all time. A similar situation can take place in the 
circuit DCE experiments, where the artificial detecting 
"atoms" do not move at all (unless some scheme of turn- 
ing on their interaction with the field at precisely chosen 
instants of time is used). 

Although this second scheme can result in diminish- 
ing the number of created photons flEl . [l6l | , it can have 
some advantages from the point of view of the photon de- 
tection and generation of novel quantum states, as was 
shown recently in [13, HH . But in that papers the prob- 
lem was treated only numerically, since the exact Rabi 
Hamiltonian does not allow for simple analytical solu- 
tions. At the same time, it is desirable to have also at 
least approximate analytical solutions, which could help 
us to understand better the mechanism and details of 
the process. Finding such solutions and comparing them 
with numerical ones is the goal of this report. 



(in some wide sense). Then the number of created pho- 
tons can be, in principle, even bigger than in the cavity 
DCE (since the time-dependent part of the Hamiltonian 
becomes linear with respect to the annihilation/creation 
operators, instead of the quadratic Hamiltonian ([l}, as 
was mentioned in [3]). However, namely the cavity DCE 
seems the most impressive, from our point of view, and 
for this reason we concentrate on it. Formal mathemat- 
ical relations between the cases \uj(t), g = const] and 
{w = const, g{t)} were discussed in [20|, where the influ- 
ence of dissipation was also taken into account. 

The first step in obtaining analytical solutions is to go 
to the interaction picture by means of the time-dependent 
unitary transformation 

\^{t)) = Vitmt)), V it) = exp [-It (77/2) in + a,/2)] . 

Then the interaction Hamiltonian acting upon the new 
function \ilj{t)) reads 



H = V^t)HoV{t)- iV^t)Vit) 

= giaa+ + aV_) + iqia)"^ — a?') — xn 



(4) 



II. ANALYTICAL SOLUTIONS FOR 
SIMPLIFIED HAMILTONIANS 

Following the studies d, we choose the time de- 
pendence of the cavity frequency in the harmonic form 
bjit) = a;o[l + esin(77t)], where e ^ 1 is the modula- 
tion amplitude and 77 is the frequency of modulation. It 
is known J^] that in the absence of the atom-field in- 
teraction the mean number of photons grows with time 
exponentially if 77 w 2u)q and eoj^t > 1. Hereafter we 
normalize the unperturbed cavity frequency to 1, writ- 
ing the modulation frequency as 77 = 2 (1 + x), where x 
is a small resonance shift. Moreover, for e ^ 1 we write 
cj (i) ~ 1, as the modulation influence is only relevant for 
the squeezing coefficient X (i) ~ (£77/4) cos(77t). Having in 
mind obtaining approximate analytical solutions, we re- 
place the Rabi Hamiltonian ^ by its Jaynes-Cummings 
reduced form f23| Hjc = 5 (a(T+ -f a V_ ) . Therefore, 
our starting point is the Hamiltonian 

Hq =n + ixit)ia)'^ ~a'^) + ^(j^+giaa+ + a)(j^) (3) 

in the weak atom-field coupling regime, \g\ <^ l,ri. We 
wish to find approximate analytical solutions to the time- 
dependent Schrodinger equation |vE'(t)) = — ii?o|^(^))- 
This problem was solved in [15[ under the restriction 
|e| <C \g\, but now we consider the case of arbitrary (al- 
though small) values of e and g. 

Note that we assume that x(0 the only time- 
dependent function, while the coupling constant g is 
time-independent. In principle, photons can be created 
also in the case of fast variations of the coupling con- 
stant g, instead of the cavity frequency uj 23]. This case 
can be also interpreted as another realization of the DCE 



where q = £{1 + x) /A and A = 1 — is the detuning 
parameter. In the following we consider the separable 
initial state with the field in the vacuum state, and the 
atom in an arbitrary pure state 



IV'(0)) = (a|.g)+/3|e))®|0) 



/3 = ^/l-a^, 



(5) 



and for simplicity we assume real a and /3. Explicit (al- 
though approximate) analytical expressions for the prob- 
ability amplitudes and the average values of the main 
system observables can be obtained in two regimes: the 
dispersive and resonant ones. 



A. Dispersive regime 

The dispersive regime occurs when |A| ^ g. In this 
case we can simplify Hamiltonian ^ by means of the 
unitary transformation H' = U HU'^ with 

C/ = e^^, Y^a)a^-aa+, C = 5/A < 1. (6) 
Expanding exp (C^^) in the Taylor series, we obtain 
H' = iO (a)"^ — a^) — ipn — thOz — 2iqC (a^cr+ — a(T_) 
-^-l9C^cina++na^a-) +gO{C'') , 
where 6 = <C 1 is the dispersive shift and 

ip = x + 5c7^, 6* = q(l + C^cr^), w = iA + S + x)/2 . 

The DCE is described by the term i9 (a'^^ — a^). In the 
absence of other terms, it would result in a slow evolu- 
tion of the photon operators (in the Heisenberg picture) 
on the time scale of the order of e~^. On the other hand. 



3 



the term wdz alone would result in oscillations of the op- 
erators a± with the frequency 2w. Consequently, under 
the condition |£| ^ |2t(7| the terms containing operators 
CT± can be considered as rapidly oscillating (and small 
due to the presence of coefficients C or ) , so that they 
can be removed following the standard ideology of the 
rotating wave approximation. In this way we arrive at 
the following effective Hamiltonian, which takes into ac- 
count the terms up to the second order in ^ (we neglect 
the unessential constant term —5/2): 



Lpn — ZUCTz 



(7) 



It is valid roughly for \5\t <C 1. Then the time-dependent 
wavefunction (in the interaction picture) can be written 
as IV'(O) — U'^Kip,eeyi^{ivjt(Jz)U\'ip{Q)), where we define 
the squeezing operator 



A,, 



exp { \ivin + V2 (a^^ — a^)] 



(8) 



Note that ip,6,vi,V2 are operators with the respect to 
the atomic basis, containing the diagonal matrix az- The 
operator A^^^vi has the following properties [2^ : 



-ivit/2 / c 



2"n! 



where \k) denotes the k-th Fock state and 
Cvx,v2 = cosh(d„^^„2t) - idj^^^^^wi sinh (d„^^„2t) , 



Svi,v2 = 2d„j^„2W2 sinh(d„i^^,2i) , d^^^^^ = y^4u| - w^. 



For ui = we use the shorthand notation aJ)"-* = Aq"^^^, 

= Co, = cosh(2u2t), = Sq^v2 — sinh(2w2i)- 
After some manipulations one can obtain the probability 
amplitudes (exact to the second order in C) 



{e,2n\m) 
(g,2n+l|V(i)) 



(e,2n- l|V'(i)) = ae-*^*Cv2nA 



A 



+e--*C' (2n + l)C;_\,_Aj,"_) , 



,(») 

ip-,0- 



where 9± — q{l ± C^) and if± — x ± 5. Other relevant 
quantities as functions of time are as follows: 



(a) = 



-2C^( 



1 - Re 



(Eo + SO/C, 



} 



p. 



-2C2(l-Re[(I]o + I]i)/C;_^,_])}, 



where ^ e^-* ^^=0 (2^)' A^-k ^S.^^ ' 

We see that the exponential photon generation from 
vacuum is possible only if > (accordingly with 

the atomic initial state), therefore the resonance shift x 
must be adjusted as function of the atomic initial state. 
Thus, without the resonance adjustment, i.e. setting x — 
0, there will be no exponential photon creation for 9± < 
\S\ /2, or roughly for £ < 2 \S\. Moreover, ioi x = 6 (or 
X — —5), only the initial state a\g, 0) (/3|e, 0)) contributes 
to the exponential photon growth whenever e < 4 

In particular, for the initial state |(7,0), by adjusting 
the resonance shift to the most favorable condition x = S, 
one obtains the following simple expressions: 



(n) = (1 - C^) smh^{2e^t), 



(9) 



Q = l + 2(7^)-C^ 5^2) = 3-|-sinh-2(26'_i), (10) 



2\ „46l_t 



(11) 



((Am = :^[(l-C 



2\ „-49_t 



+ C'] , (12) 



where Q = [((An)^) — {n)]/ {n) is the Mandel factor 

and 5'^^ = (fltataa) / (rt)^ — 1 + Q/ {n) is the second 
order coherence function. The field quadratures (in the 
interaction picture) are defined as 

X=(a + a^)/V2, P= (a-at)/(^/2i). 



B. Resonant regime 

In the resonant regime, A = 0, we choose the unitary 
operator U in the form 



U = e'^^, Y = a^a++aa_, ^ = 2g/e, (13) 

assuming the weak coupling regime, ^ ^ 1 (in the op- 
posite case, ^ 3> 1, no more than two photons can be 
created due to the fast exchange of excitations between 
the field and atom [15]). As shown in [lO, HH, in this 
case it is convenient to set x = 0, so one obtains 



35,. 



H' 



{ana+ - a)'^a+ + h.c.) + 0{t 
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where h.c. denotes the hermitian conjugate. Thus, to the 
second order in ^ we obtain a shghtly modified squeezing 
effective Hamiltonian 

vahd roughly for gt ^ 1. Then the wavefunction is 

|V(t)> = C/tAo,^[/|V^(0)), 

so that the probability amplitudes to the second order in 
£, read (here d± = (e/4)(l ± £_^) and t}o = 



1) A, 



+?2(2n + l)C-;A("j} 



(g,2n-l|VW) 
(5,27.1^-^) 

{e,2n\m) 
(e,2n + l|V(t)) 



The expressions for the average values of the main ob- 
servables are as follows: 



{X) - 0, (P) - \/2C«/3 [6-2"+* _ C-^'^ 



we can believe that they are correct on the time scales 
t ^ in the dispersive regime and for t <C 1/g in the 
resonance case (so that in both the cases the product et 
can be bigger than unity, thus enabling the generation of 
many photons). To check the validity of analytical for- 
mulas obtained, we solved the Schrodinger equation nu- 
merically for the given initial conditions, using instead of 
the approximate Hamiltonian ^ the exact initial Hamil- 
tonian Hoo + Ha+Hji [i.e., taking into the account the in- 
teraction in the complete Rabi form The numerical 
results turn out to be in a very good agreement with the 
analytical ones. For example, in Fig. [l]we show the ex- 
act dynamics for the initial state ([5]) with /J^ = 0.3 in the 
resonant regime with parameters A = x = 0, g — 5-10^^, 
and £ = 2-10~^. A small noticeable difference can be seen 
only for the evolution of the probability Pe{t) of finding 
the atom in the excited state for st > 3 (the variations of 
this quantity are small due to the smallness of parame- 
ter ^ = 0.05). The analytical results for other quantities 
are indistinguishable from the numerical ones within the 
thickness of lines used in the plots. The mean number of 
photons grows exponentially for et > 2, and the photon 
statistics is super-Poissonian, since the quantum state is 
close to the squeezed vacuum state (small nonzero prob- 
abilites of the odd numbers of quanta arise just due to 
the atom-field interaction). 



I S2^_ + C 



-3/2 



2 



7-1/2 



a'[S^,_+2e 



Pe = 2a^^^ 



1 + 



-3/2 



n p 3/2 



In this case simplified expressions can be obtained for the 
initial excited state (a = 0): — \ — ^ {n), while (n), 
0, g'^^\ {{AXf) and {{APf) are given by Eqs. 
with substitutions C ? and 0^ — ?> t5+. 



III. DISCUSSION AND CONCLUSIONS 

Since the results of the two preceding sections were de- 
rived in the frameworks of small-parameter expansions. 




FIG. 1: Exact numerical results in the resonant regime and 
the analytical results for Pg. The photon number distribution 
Pn = \ {n\ip)\'^ is evaluated for et = 3.9. 

In Fig. [2]we show analogous results for the initial state 
([5]) with = 1/2 in the dispersive regime with parame- 
ters X = 0, g = 5 ■ 10~^, A = 12g, and s = 3S, comparing 
numerical and analytical values for the photon number 
distribution. Again, the coincidence is very good, since 
the differences are seen only for very low probabilities, 



5 



less than 10 ^. For other quantities, analytical results 
are indistinguishable from numerical ones in all plots. 
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FIG. 2: Exact numerical results in the dispersive regime. P„ 
is evaluated for et = 2, both analytically and numerically. 

In conclusion, we obtained closed analytical expres- 
sions for the atom-field dynamics generated by the dy- 
namical Casimir effect in the resonant and dispersive 
regimes, for arbitrary pure atomic initial state. Our re- 
sults are exact to second order in the small parameters 
C = g/A and ^ — 2g/e, being in good agreement with 
numerical data, so they can be used to quantify the in- 
fluence of the atom on the DCE for different modulation 
frequencies and atom-cavity detunings. An interesting 
unsolved problem is to extend the time interval of valid- 
ity of results from the scales of the order of or 
up to the scales of the order of or (in the ab- 
sence of field-atom coupling some results were obtained 
in ^M)- 
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